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Abstract 

The generality of the Bethe's two dipole model for light diffraction through a subwavelength 
aperture in a conducting plane is studied in the radiation zone for coated conical fiber tips as those 
used in near field scanning optical microscopy. In order to describe the angular radiated power 
of the tip theoretically, we present a simple, analytical model for small apertures (radius < 40 
nm) based on a multipole expansion. Our model is able to reproduce the available experimental 
results. It proves relatively insensitive to cone angle and aperture radius and contains, as a first 
approximation, the empirical two-dipole model proposed earlier. 
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Introduction and experimental motivation - The well-established Fresnel and perturbation 
methods used in optics are not adapted to diffraction phenomena implying strong near-field 
patterns which occur for wavelengths A comparable to the dimensions of the diffracting 
object. This is due to the fact that for small objects, the "interaction zone", of typical 
dimension A, invades the space around the diffracting matter. As a consequence, the in- 
teraction of light with matter cannot be considered as a small perturbation. Rayleigh[l|, 
Bethe[2|], and Bouwkampjs), 0] have studied the particular case of diffraction by a small 
aperture in a perfect metallic screen. This problem - hereafter referred to as the Bethe's 
problem - is very typical because it contains dipolar terms which dominate in the "radia- 
tion" far field of the hole. The present work addresses the validity and the survival of the 
diffraction phenomenology revealed by the Bethe's model in the case of small aperture at 
the apex of a conical screen. 

Such an unusual geometry is encountered in Near Field Scanning Optical Microscopy 
(NSOM). Indeed in NSOM, a conical optical tip with a sub- wavelength aperture at the 
apex is raster-scanned in close proximity to the investigated surface. Most NSOM make 
use of the tapered and coated optical fiber introduced some 10 years ago by Betzig et al. 
[5l. The angular radiation (far field) profile of the tip depends on the aperture diameter as 

n q 

has been extensively discussed by C. Obermuller and K. Karrai |6J, |7|. For small aperture 
sizes, the far-field radiation pattern is polarization-dependent with a large backward emis- 
sion in the P-polarization in contrast to the S-polarization (incident electric field parallel and 
perpendicular to the analysis plane, respectively) for which there is little backward emis- 
sion. For small radiation angles, a simple combination of mutually perpendicular electric 
and magnetic radiation dipoles situated in the aperture plane reproduces the main features 



of the angular radiation pattern for the smallest apertures (60 to 80 nm diameter) 
There is presently no formal justification for this 2-dipole model, and for that reason, the 
transmission profile of tapered optical probes, i. e. the problem of light diffraction through 
a subwavelenght hole in a conical screen, needs further theoretical examination. 
In this letter, we solve Maxwell's equations analytically in the radiation zone with the use 
of appropriate boundary conditions in order to get insight into the peculiarities of a conical 
diffraction geometry compared to the more common planar geometry. As a critical test of 
our analytical model, we compare calculated emission profiles with available experimental 
profiles for both polarization states. Our analysis, based on a conical tip geometry with an 
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aperture radius that is negligible compared to the wavelength of the incident light, is able 
to reproduce the radiation pattern of small apertures over a large angular range. Remaining 
discrepancies at the largest radiation angles can be accounted for by the boundary condi- 
tions in the metallic coating. The above-mentioned empirical radiation dipole approach is 
formally justified for small radiation angles. In addition, it is shown that the angular power 
spectrum does not critically depend on the exact tip geometry. Especially the dependence 
on the taper angle is found to be very weak. 

Scatterring and general properties - In order to illustrate the main physical effects of small 
apertures in metallic surfaces, it is worthwhile to look at similar problems that have been 



treated historically. As already mentionned in the introduction, Rayleigh 1| and Bethe 
(see also Ref. [3]) have studied the light diffraction by a subwavelength aperture in infinite, 
perfectly conducting screen. They have demonstrated that, in the far field, the diffracted 
electromagnetic field can be described by a combination of effective dipoles which depend 
on the local incident field (Ei nc ,B; nc ) and on the hole size a: 

a 3 3 
Pe ff =-E inc , ± , M eff =-2^B mCi|| . (1) 

The symbols _L and || refer to the perpendicular and tangential components of fields, respec- 
tively. These dipole terms are very similar to the induced dipole terms which appear in the 
well established Rayleigh-Mie theory of diffraction by a small perfectly conducting sphere: 
Pinduced = a 3 Ei nc , M in d U ccd = — ^rB inc . The principal difference arises in the substitution 
2 — > ~ between the magnetic and electric dipole. Strong analogies exist between the theory 
of diffraction by small hole and the theory of diffraction by small particles. We can under- 
stand the physical origin of these similitudes using the Clausius Mossoti formula. Indeed, 
following these formula, the polarization produced by a locally constant electromagnetic 
field E , B in a small dielectric sphere of constant permittivities e, \x which is immersed 
in a homogenous medium with permittivities eo,/io is defined by: P = eo -7~^ Q 3 Eo, and 
M = ~ ^°^2 fl3 ^ 0. With the condition e/eo, /-to//-* <C 1 for a hole in a perfect metal, we 
obtain the following two relations: P = — ^-a 3 E , M = ^a 3 B which are related by the 
same factor of 2 as for the Bethe case. In the opposite case, namely if e/eo,/io/V ^> 1, we 
obtain the Rayleigh-Mie's solution. 

The diffraction of light by a small aperture in a perfectly conducting plane can be qualita- 
tively regarded as the complementary case to the scattering of a wave by a small conducting 
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particle. We can then consider that the incident light "sees" the hole as a (P,M) combi- 
nation. The general form of Eq. [TJ is certainly conserved for various quasi circular opening 
in planar geometry and must essentially be independent on the aperture shape of subwave- 
length holes. As a good staring point we can anticipate the generality of this result, and we 
can admit that in any problem of diffraction by subwavelength holes in non planar screens 
the emission zone can be described by dipoles terms similar to Eq. [TJ Nevertheless, numerous 
differences certainly arise from boundary effects when we look at the field near the metallic 
surface. In particular, the dipole orientations may change with the screen geometry as we 
will show later on. 

Model - In the case of an NSOM fiber tip, the aperture is situated on the apex of a conical, 
metal-coated structure (see Fig. [TJ. For eigenvalue-problems of conical symmetry around 
some origin, it is convenient to express the fundamental solutions in spherical coordinates 
(r,8,(p). The separation of the angular and radial variables in the Helmholtz equation fol- 
lows the "Hankel(r)-spherical-harmonics(#, 0)" solutionf? ] ip (r) = h u (kr) ■ Y UjTn (8, 0) with 
the modification that in this geometry the parameter v can have a non-integer value which 
is determined by the boundary conditions. In analogy to the method of multipole expan- 
sions with spherical harmonics used in free space, we can develop our propagating field as 
an expansion in electric and magnetic "quasi-multipoles" . The principal difference to the 
classical multipole expansion {9] consists in the substitution I — > v. 

We find that in general, a propagating electromagnetic field in conical geometry produces 
an angular distribution of radiated power of the form 

+ H)- +1 ]T aZ] m ■ LF -; m(M) || 2 , (2) 

with L = \t x V. Here, z/g and vm refer to the electric (TM) and magnetic (TE) quasi- 
multipoles, respectively. The coefficients a UtTn depend on charge and current distribution on 
the metallic cone. It can be shown that for two roots v and v' (with v < 1/'), the ratio a v i ja v 
is comparable to {ajX) v ~ v 1 which can be neglected. Therefore, we can omit all terms 
in the expansion except for the first two roots ve,o, and vm,o- 

Specific boundary conditions can be written for a good metal such as aluminum. They 
depend only on its electric permittivity e (u) (~ —34.5 + 8.5i for A = 488 nm). If we 



neglect the skin depth 5 (u>) we arrive at the following boundary conditions on the cone: 
9Y UMtin (9, </>) /dcos9\s — 0, and Y VE ,m (9, 4>) \s — 0. However, we can make a more realistic 
approximation and use a finite value for the skin depth in the optical domain, (5 ~ 10 nm). 
In this case, there exist just outside the surface tangential Eii and By linked by the relation 
E|| = a/ /u/eB|| x n. In this equation n is inwardly directed normal on the metal surface and 
yjife can be interpreted as a surface impedance. This relation determines an eigenvalue 
condition on v and authorizes an infinite but discrete number of solutions which depend 
only of the tip angle /3. For a value of (5 on the order of 10°-15°, we obtain ue,o — 0.95 
and i>Mfl — 0.96. The v values and, consequently, the angular profile are not sensitive to 
variations of the tip angle unless it exceeds 20°. 

It is important to note that the boundary conditions cannot be applied inside a sphere of 
radius a/ (sin j3) because it contains no metal. This so-called aperture zone represents the 
"terra incognita" near-field zone (r/A < 1) of the fiber tip (cf. Fig. [1]). The coefficients 
a^ m can be approximated by a Taylor series around the tip origin, m = ±1 are the only 
authorized terms for linearly polarized light propagating down the fiber. These terms possess 
the form a E ~ P e // • E , a M ~ M e // • B which depend on the field E , B = z x E at the 
aperture and on the two dipoles: 

1 a 3 



eff 



r E , M eff = 2z x P eff . (3) 



4 a/2 sin /3 3 

The normalized angular distribution of the radiated power can then be expressed as 



1 dP „ 2 
Pdfi ~ 5 



Tm = == x e r + 2Re 



(4) 



v E \Ve + 1) \Jv M \vm + 1 

In this formula Re and Im are related to the real and imaginary parts of the generalized 
harmonic functions. It appears also a characteristic factor 2 between the electric and the 
magnetic quasi-multipole. It is important to remark that neither the absolute value of the 
radiated power nor the absolute value of the dipoles are described by our model. This would 
involve a much more elaborate description of the core/cladding characteristics inside the 
fiber tip. We can also note that the normalized profile does not depend on the aperture 
radius. Long wavelength light is blind to subwavelength details, and therefore the aperture 
size does not constitute a critical parameter of our model. 

Eq. H] reduces to the Obermuller-Karrai model when we put v — 1 (a value close to the real 
values ve,o and z/ M0 )- This corresponds to the ideal case of a linear, radiating two-dipole 



antenna. In this limit f3 and a//3 tend to zero, and we can write 



1 dP v - 1 _ 3 || (e r x P eff ) x e r + M cS x e r 
Pdn = 8^ ||P cff || 2 + ||M cff || 2 




The effective electric and magnetic dipoles are linked by M e g = 2e z x P eff and by the 
generalized Bethe relation P eff oc a 3 E inc where E inc denotes the incident light polarization. 
Also, the factor of 2 reappears which is common to all cases of diffraction by a subwavelength 
aperture as precedently explained. The two-dipole model has another significance because 
the general expression Eq. 0] converges to Eq. [5] for small azimuthal angles 9. As a direct 
consequence the phenomenological Obermuller-Karrai's model is justified and constitute a 
good approximation valid for small azimuthal angles. 

Results and discussion - The following Figs. 121131 show experimental data for both S- and 
P-polarization and compare them with the calculated profiles of our model. Also shown are 
the curves for the simpler two-dipole model which is a good approximation only for smaller 
angles. There are no adjustable parameters in both models. 

As shown in Fig. [2] (S-polarization), our model and the Obermiiller-Karrai model are 
equivalent for the experimentally accessible angles. They follow the experimental data for 
all recorded angles. This behavior is different for the P-polarization as shown in Fig. |3j 
Our model predicts an important "backscattering" in good agreement with the experiment. 
The two-dipole model does not correctly predict this effect. Qualitatively, the difference 
between both polarisations can be understood by realising that the Poynting lines in the 
P-polarisation are tangential to the metal surface, and are consequently very sensitive to 
boundary conditions on this surface, while in the other polarisation, Poynting lines tend 
to avoid the metal surface because both the electric and magnetic fields vanish there. The 
quasi-multipole model is nevertheless in disagreement with the measured data for extreme 
9 values. We believe that this discrepancy is due to several effects. First, there are strong 
near-field effects on the surface of the metallic cone which are not accounted for in the 
model. They appear in higher order terms in the quasi-multipole expansion (these terms 
are not considered here). Additionally, incoherent scattering due to surface roughness 
reduces the collected light for large angles which could explain the discrepancy. 



Conclusion - In summary, we have developed an analytical model for the far-field trans- 
mission pattern of conical, metal-coated fiber tips with a small, subwavelength aperture. It 
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FIG. 1: Geometry of a NSOM fiber tip. It is characterized by i) the aperture radius a in the metallic 
coating and ii) the tip angle f3. The quasi-multipole model is valid only in the far-field zone where 
a/(sin/3) < r. The boundary conditions cannot be applied in the near- field (or aperture) zone 
because it contains no metal. 

contains, as a first approximation, the two-dipole model anticipated earlier [6]. It justifies the 
generality of the Bethe's result for the case of a conical screen. Our model is in good agree- 
ment with the experimental data and explains in particular the backward emission found 
for the P-polarization. The angular transmission profile does not depend on the aperture 
radius for small holes, and is insensitive to the tip angle as long as it is smaller than 20°. 
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FIG. 2: Normalized angular distribution of radiated power in S-polarization where the detector 
is scanned perpendicular to the plane of polarization of the incident light. Experimental data are 
shown in circles. The thick line corresponds to the present quasi-multipole model, the thin line to 
the 2-dipole model. The experimental data (aperture radius a = 30 nm, light wavelength A = 633 
nm) are taken from Refs. 2,3. 
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FIG. 3: Same as Fig. [2] but for the P-polarization where the incident electric field is parallel to the 
analysis plane. 
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